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Abstract. We present the analytical solution of a maser 
shaped like a spherical shell. We determine the general 
condition on the size of the central cavity at which a sphere 
becomes a shell maser, and derive the intensity, beaming 
angle and observed size of both unsaturated and saturated 
shells. 
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1. Introduction 

Given a pumping scheme, the brightness temperature of 
an unsaturated maser can be immediately determined 
from the source's length along the line of sight because 
the gain, r, is simply proportional to it; there is no need 
to even consider the geometrical shape. In a saturated 
maser, on the other hand, the dependence of gain on dis- 
tance traveled through the source varies between the satu- 
rated and unsaturated domains. As a result, the brightness 
temperature in any given direction cannot be determined 
before a complete solution of the full maser structure has 
been attained. An analytic solution of this highly non- 
linear problem is possible because, in all sources with an 
appreciable gain, maser radiation is highly beamed (e.g. 
Elitzur 1990b; hereafter E90): at any location in an ar- 
bitrarily shaped maser, the longest chord through the 
source provides the direction of the local dominant ray 
and maser radiation is mostly confined to a small beam 
around it. The complete, detailed solution is highly de- 
pendent on the source geometry; indeed, it cannot even 
be attempted without fully specifying the geometry. Thus 
far, detailed solutions have been derived for linear masers 
(Elitzur 1990a), for the planar geometry of a disk (Elitzur 
et al. 1992), and for the three dimensional configurations 
of spheres and filaments (Goldreich & Keeley 1972, E90, 
Elitzur et al. 1991). 
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Fig. 1. Geometrical notations for a shell maser. 



Recent interferometric observations of masers in late- 
type stars reveal ring-like structures, indicative of a shell 
geometry. The observations include both H2O (Reid & 
Menten 1990) and SiO (Diamond et al. 1994, Miyoshi et 
al. 1994, Greenhill et al. 1995). The shell geometry is fun- 
damentally different from all the maser geometries con- 
sidered until now, which will be referred to as "regular" . 
In regular geometries there is a single dominant ray at 
virtually every point (the only exception is an occasional 
singularity, such as the center of a sphere or a disk maser). 
The dominant ray is aligned with the local axis of sym- 
metry and its direction coincides with that of the flux 
vector. In contrast, at an arbitrary point in a shell there 
are many chords of maximal length belonging to the rays 
that graze the core (see Fig. 1). Thus the dominant rays 
in a shell maser span the surface of a cone whose axis is 
aligned with the radius vector. From symmetry, the flux 
vector points in the direction of that axis and it does not 
coincide with the direction of any dominant ray. 

Maser emission from shells has been studied thus far 
only in numerical calculations (Western & Watson 1983; 
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Bujarrabal 1994). Utilizing the techniques developed in 
the studies of regular geometries, we present here the com- 
plete analytical solution of shell masers. Since the solution 
methods are described in detail in the referenced litera- 
ture, the presentation here centers mostly on the features 
unique to the shell geometry. 

2. Generalities 

Consider a maser shaped as a spherical shell with inner 
radius Ri and outer radius R (Fig. 1). As in all previous 
solutions, the shell is assumed uniform and quiescent (no 
ordered motions). At a point in the shell characterized by 
radius r, the longest chords belong to the rays that graze 
the core, inclined at 9o = sin -1 (Ri/r) to the radius vector. 
Introduce dimensionless variables x and Xi such that 



r = xR . 



Ri — XiR , 



#o = sin 1 (x i /x). 



That is, Xi is the relative size of the cavity, and positions 
in the shell are characterized by Xi < x < 1. 

Denote the intensity along a dominant ray of any maser 
la. In regular geometries the radiation is confined mostly 
to 1 9 1 < 8b, where 9 is measured from the local symme- 
try axis and 9b (<C 1) is the beaming angle. Therefore, 
the angle-averaged intensity is J ~ \lo^b' ^ the inten- 
sity is assumed approximately constant inside the maser 
beam. In shell masers, on the other hand, the radiation is 
confined mostly to 9 < 9 < + 8b (see Sec. 3 below). Ne- 
glecting rays outside the main beam, the angle-averaged 
intensity for a constant intensity inside the beam is 



Isin9d9 ~ \l®9b (sin0 o + \®b cos O ) • 



(2) 



Instead of a rectangular intensity profile we can employ 
an exponential fall-off, / oc exp[— (9 — 0o)/0&)]- This angu- 
lar distribution produces the same result except that the 
factor | inside the parentheses is replaced by unity. 

When 9 <C \9b (<C 1), the second term in the paren- 
theses dominates and the result of regular geometries is 
recovered. Therefore, the shell is properly described by 
the spherical maser solution as long as its central cavity is 
sufficiently small that 9a < \8b- The shell structure must 
be considered only in the reverse situation in which the 
first term in the parentheses dominates and the spheri- 
cal solution no longer applies. So the transition to a shell 
maser occurs when 



9q > 9b , 



(3) 



omitting the factor i to accommodate also the exponential 
profile. This condition, similar to the filamentary condi- 
tion that defines elongated masers (Elitzur et al. 1991), 
will be referred to as the shell condition. At the inner 
surface (x = x{), 9 = n/2 and the condition is trivially 
obeyed. It provides the most meaningful constraint at the 
outer surface (x = 1), where 9q is minimal. 



The meaning of the shell condition is quite simple. On 
the surface of a spherical maser, the radiation fills a cone 
with an opening angle 9b centered on the radius vector. 
Now begin to carve a hole at the center of the sphere. As 
long as that hole is contained in the radiation cone it can 
be ignored and the maser remains a sphere. The maser 
becomes a shell once the hole emerges from the side of the 
radiation cone, i.e., when the shell condition is obeyed. 
The angle-averaged intensity in a shell maser is thus 



J = ^I 9bsin9 



llo9 b — ■ 

1 x 



(4) 



To a distant observer, the maser appears as a ring-like 
structure with inner radius Ri and outer radius Ri + H, 
where 



(1) H = RB b Jl-x 



(5) 



and 9b is the beaming angle at the surface. Observations 
can directly determine only the inner radius Ri and the 
thickness H. The outer radius R is unobservable, thus the 
ratio Xi cannot be directly determined. This problem ex- 
ists of course in all regular geometries, too, where the only 
maser dimension amenable to observations is the observed 
size in the plane of the sky, not the length along the line 
of sight. Neither the radius of a spherical maser nor the 
length of a filamentary maser nor the length along the 
line of sight of a planar maser are directly observable (e.g. 
Elitzur ct al. 1992). 

Consider now the maser flux. In regular geometries 
the flux vector points in the direction of the dominant 
ray and is dominated by the contribution of 9 = 0, thus 
its magnitude obeys F ~ Air J. In a shell maser, on the 
other hand, the directions of the dominant ray and flux 
vector are different. Since the flux is dominated by the 
contribution of 6 = 6o, its magnitude obeys 

F = J cos Oldfl ~ cos O J Idfl 

= cos6> x47rJ = 4nJy/l - (xi/x) 2 . (6) 

The result for regular geometries is recovered when Xi = 
0. The flux vanishes on the shell's inner surface because 
non-dominant rays were neglected. 

3. Unsaturated shell masers 

The optical depth in an unsaturated maser is always pro- 
portional to pathlcngth. At an arbitrary position in the 
shell, denote the length of the longest chord £q. A ray in- 
clined to this chord at a small angle 59 (> 0) has a length 



r cos(6> + 59) + \J R 2 — r 2 sin 2 (0 o + <50) 
t £ - £i59 - 



e 2 59 2 , 



(7) 
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where 
4> = 

li = 

4 = 




(8) 



In regular geometries the dominant ray corresponds to a 
true maximum of the geometrical surface, i.e., d£/d9 = 
and the variation of £ with 59 is quadratic. In a shell the 
dominant ray does not correspond to a true maximum, 
the first derivative does not vanish and the variation of £ 
with 59 includes a linear term (which properly vanishes 
when Ri = 0). As a result, the angular distribution in an 
unsaturated shell is 

1(69) = Ioe -KohS9-Ko£ 2 S9\ (9) 

Here kq is the unsaturated absorption coefficientB and 
lo = Sexp(noio), where S is the source function. The 
intensity detected by an external observer is obtained 
by substituting x = 1, which produces £o = 2£ 2 = 
2R^/l — x\ and £\ — 2Ri. Instead of the standard Gaus- 
sian that typifies the angular distribution of spherical 
maser radiation, the argument of the exponential contains 
a term linear in 59. At very small 59 the linear term dom- 
inates, at very large 59 the quadratic term. The shell con- 
dition ensures that the linear term dominates inside the 
entire maser beam, 59 < 9b- With the dominance of this 
term, the radiation at any point in an unsaturated shell 
maser varies with direction in proportion to exp(— 59/9b), 
where 



9,, 



1 



K £l 



(10) 



is the beaming angle. At the shell's outer surface 9b = 
l/(2KoRi), and the shell condition becomes 



1 



kqR 



(11) 



A factor of \ was omitted on the right-hand side to provide 
a slightly more stringent constraint, in line with another 
estimate described below. 

Only rays with 59 > were considered thus far. Rays 
with 59 < cut through the cavity and the decrease in 
their length from f is dominated by a term proportional 
to \/—59, much steeper than the linear decline for 59 > 0. 
For all practical purposes, the angular distribution can 
always be assumed to have a sharp cutoff so that only 
9 > 9q need be considered. 



1 Following what is by now a standard procedure, the nega- 
tive of the absorption coefficient is used for the maser transition 
so that it is positive. 



The condition just derived was determined for unsat- 
urated shells. However, a maser cannot become saturated 
unless it was first unsaturated and had to obey this condi- 
tion during its unsaturated phase to be considered a shell. 
Therefore, Eq. [ll] is the general condition that must be 
obeyed by all shell masers. For saturated shells we have 
to consider the magnitude of k R that the maser must 
have had to become saturated. Typically, saturation oc- 
curs at KqR <^ 10, thus astronomical masers obey the shell 
condition when Xi <^ 0.3. 

General results for arbitrary Xi, displaying explicitly 
the transition from a sphere to a shell, can be obtained 
from the full angular distribution of Eq. ^. The require- 
ment that the intensity drop to 1/e of its peak value pro- 
duces a quadratic equation for the beaming angle whose 
solution is 



- 1 



(12) 



where z — kq£i / {2\J k £ 2 ) ■ When z 3> 1, equivalent to 
the shell condition of Eq. O, the beaming angle of Eq. [l(] 
is recovered. The opposite limit, z — » 0, properly re- 
produces the beaming angle of an unsaturated sphere, 
9b = 1/y/noR. Furthermore, the complete expression for 
J can be calculated from an integration of the full angular 
distribution, and the result is 



1 



J = jlo — 7T 1 _ n/tt(z - 9q\/ n Q £ 2 )e z erfc z 
kq£ 2 L 



(13) 



Again, utilizing properties of the error function erfc at 
large and small z reproduces the proper expressions for J 
in the respective limits. 

From Eq. |^, the observed thickness of the maser ring 
obeys 



k H = 



(14) 



This result can also be obtained by noting that the surface 
brightness observed at infinity varies as e~ h / H , where h is 
the parallel linear displacement from the brightest ray. It 
is interesting that the observed thickness of unsaturated 
shells is independent of size,cl the shell dimensions enter 
only in the ratio Ri/ R. 

Since the thickness H was derived assuming the shell 
condition, it diverges when xi — > 0. Employing the full 
expression for the beaming angle, Eq. PL rem oves the di- 
vergence and properly produces H = \J R/ kq when xi = 
0, the observed size of an unsaturated spherical maser. 
In addition, the area generally provides a useful indicator 
for the effects of geometry on dimensions of the observed 
maser (see table 1 in Elitzur et al. 1992). The observed 



2 The maser size enters only indirectly through the condition 
Kofi > 1, required to ensure that the maser provides amplifi- 
cation sufficient to cause beaming. 
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area of an unsaturated shell maser, where A b s — 2irRiH, 
is 



-4 



obs 



2k q ' 



(15) 



where Iq is the longest chord at the outer surface. By 
comparison, in an unsaturated sphere, where £q — 2R and 
the observed region is a small cap, A b s = 2tt£q/kq. Apart 
from a factor of 4, the expressions for the observed areas 
are the same. 

4. Saturated shell masers 

Consider now a sequence of masers with a fixed a;, , increas- 
ing n n R. If the original maser obeys the shell condition, 
Eq. |ll], subsequent masers obey it with a larger margin. 
The emerging intensity grows exponentially with maser 
size and at a certain radius, R s , it obeys J(x = 1) = J s , 
i.e., the shell saturates. The corresponding radius obeys 



— — exp ( 2k Ri 
k R s V 



4 7 , 



(16) 



where 7 = J s /S; typically, astronomical masers have 7 <• 
10 5 — 10 7 . The approximate solution of this equation is 



k R s 



1 



2^1" 



dn 



27 



VT^xf 



(17) 



a solution valid when the argument of the logarithm on the 
right-hand-side is large (e.g. E90). When Xi — 0, this result 
reproduces the saturation radius of a spherical maser. 

With further increase in radius, the shell develops a 
two-zone structure: the inner zone, Ri < r < r s , is un- 
saturated and has k ~ kq, the outer zone, r s < r < i?, 
is saturated and has k ~ kqJs/J- The boundary radius, 
r s = x s R, is still unspecified. It is R s for a maser that 
has just saturated and should shrink logarithmically as R 
increases. 



4-1- The saturated shell 

From Eq. 6, the flux divergence relation in a shell maser 
is 



dx 



X\ x 2 — xfj — kJRx 



(18) 



Since in the saturated shell k.J ~ kqJs, the solution there 
is immediate: 

x\J x 2 — xfj — x s \J x 2 . — x 2 J s — ^J s k R(x 3 — x 3 s ). (19) 

Note that when x 3> x s , J — ^J s nor, same as the result 
for a sphere (E90). 



Along a ray that grazes the core, I — Iq and di = d£o- 
Thus the radiative transfer equation past the mid-point of 
this path is 



xf dl 



1 h — — kRIq. 

x z ax 



(20) 



With the aid of Eq. [|, the transfer equations for J and Iq 
produce an equation for the beaming angle 



d\n(0 b x l /x) 2x 2 — xf 

dx x(x 2 — xf) 

The solution is straightforward, the result is 



ib = 



„2 ' 



(21) 



(22) 



where 9 S is the beaming angle at x = x s ; again, 2^ = 
reproduces the expression valid in the saturated zone of a 
spherical maser (E90). The angle 9 S can be obtained from 
the result for an unsaturated maser (Eq. |Io|) , yielding 



1 



1 



KqRXi y/l - xf + ^[xf 



K Ri 



(23) 



In the last relation we assumed that the thickness of the 
saturated zone greatly exceeds that of the unsaturated 
region, i.e., 1 — Xi S> x s ~ Xi, an assumption justified for 
strongly saturated shells. The thickness of the observed 
maser ring is 



H = 9. 



B 2 



and its area is 



A - 2?F 
K 



B 2 



(24) 



(25) 



By comparison, the observed area of a saturated sphere is 
■Kr s /Ko. 

The final quantity needed to complete the solution is 
r s , the radius of the boundary between the saturated and 
unsaturated zones. It can be obtained from the require- 
ment that, amplification by the unsaturated shell brings 
the intensity of the subordinate ray at core entrance to the 
dominant ray corresponding to J s upon exit. The detailed 
steps are described in E90 and will be omitted here. The 
final result is controlled by the parameter 



In 



I27Z? 



(« i?) 2 v / i 3 ^ 



(26) 



When the argument of the logarithm is much larger than 
unity, which is usually the case in astronomical masers, 
the boundary of the unsaturated zone obeys 



[KoRi 



(27) 
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This enables us to complete the solution for a saturated 
shell maser. The beaming angle of the emergent radiation 
is 



its intensity is 

/o = ^M"oR) 3 \fl-x 2 i (29) 
and the thickness of the observed maser ring obeys 

K ° H = v^w ■ (30) 

4-2. Complete saturation 

As kqR increases while Xi is held fixed, the parameter b de- 
creases logarithmically and r s — > Ri (see Eqs. |6]and |27| ); 
that is, the size of the unsaturated zone shrinks. In anal- 
ogy with core saturation of regular geometries, the unsatu- 
rated zone disappears when R exceeds a certain threshold 
R c and the shell is saturated throughout. A reasonably 
accurate estimate for the radius R c can be obtained by 
setting the argument of the logarithm in Eq. ^ to unity, 
yielding 

K #c^3.5 7 1/2 Zi(l-z?)- 1/4 . (31) 

Complete saturation requires dimensions that greatly ex- 
ceed typical sizes of astronomical masers, thus it is mostly 
only of theoretical interest. 

5. Discussion 

A shell becomes a sphere when xi — ► 0. Some of the ex- 
pressions we list are based on the shell condition and di- 
verge in that limit. These divergences can be avoided by 
employing the full expression for the unsaturated beam- 
ing ang le (Eq. ||) instead of its shell limit (Eq. ||). The 
results display explicitly the transition between the two 
geometries at the price of greater complexity. 

While a frequency dependence was never indicated ex- 
plicitly, it enters through the frequency variation of kq, 
the unsaturated absorption coefficient. The full absorp- 
tion coefficient, including saturation effects, used in this 
study was assumed to obey the standard form proposed 
by Goldreich & Kwan (1974). This implies that our re- 
sults are valid around line center up to frequency shifts 
of at least two Doppler widths (Elitzur 1994), covering 
practically all cases of interest. A more significant limita- 
tion is that, as with all maser solutions available to date, 
the one presented here is valid only for quiescent material. 
Complete analytic solutions for masers with large veloc- 
ity gradients do not yet exist for any geometry. While the 
impetus for this study comes from the SiO interferome- 
try results which indicate ring- like structures, these ob- 
servations also show that the maser emission comes from 



many compact spots, each one characterized by a differ- 
ent Doppler velocity (Diamond et al. 1994, Miyoshi et al. 
1994, Greenhill et al. 1995). Each feature may thus corre- 
spond to a clump, defined either by a local enhancement of 
the pump rate (reflecting a favorable combination of con- 
ditions such as density, abundance, temperature, etc.) or 
velocity coherence in a turbulent medium. In either case, 
the geometry relevant to the maser radiation is the local 
one of the individual features, perhaps best described by 
elongated structures. On the other hand, the significant 
amount of flux lost in VLBI measurements due to over- 
resolution indicates that about one half of the total maser 
emission is much more extended than the VLBI spots. It is 
also known from measurements with lower spatial resolu- 
tion that the total emission spreads over a region compa- 
rable in size to the ring seen in VLBI data. Indeed, if the 
distribution of a component of the emission smoothly fills 
up such a ring, it would very probably not be detectable 
in the high resolution experiments. Such extended masers 
would be described by the solution presented here. 

The absence of observational data on such possible ex- 
tended emission prevents further discussion at this time. 
As is the case for the spherical geometry with which it 
is closely associated, the primary importance of the shell 
maser solution may be in the insight it offers to maser 
behavior rather than direct comparison with current ob- 
servations. 
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